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Abstract 

One-loop QCD corrections to the full gauge invariant set of electroweak 
diagrams describing the charged current processes e + e~ — > u d /x P M (CC10), 
e + e~ — > u d e~ v e (CC20) and e + e~ — > u d s c (CC11) are computed. We 
compare the exact calculation with a "naive" approach to strong radiative 
corrections which has been widely used in the literature and discuss the 
phenomenological relevance of QCD corrections for NLC physics. 



1 Work supported in part by Ministero dell' Universita e della Ricerca Scientifica. 
e-mail: maina@to.infh.it, pittau@psw218.psi.ch. pizzio@to.infn.it 



Introduction 



The measurement of W mass to high precision is one of the goals of the NLC and 
will provide a stringent test of the Standard Model (SM) [|I], In fact, the mass of 
the W boson in the SM is tightly constrained and an indirect determination of Mw 
can be obtained from a global fit of all present electroweak data. The fit gives |§ 
M w = 80.338 ±0.040lg:o?lGeV where the central value corresponds to M H = 300 GeV 
and the second error reflects the change of Mw when the Higgs mass is varied between 
60 and 1000 GeV. A disagreement between the value of Mw derived from the global 
fit and the value extracted from direct measurement would represent a major failure 
of the SM. The NLC is expected to reduce the error on the measurement of W mass 
down to about 15 MeV 0], to be compared with an error, assuming design luminosity, 
of about 50 MeV at Lep 2. 

In order to extract the desired information from WW production data, theoretical 
predictions with uncertainties smaller than those which are foreseen in the experiments 
are necessary. This requires a careful study of all radiative corrections which have to 
be brought under control. In this note we will be concerned with QCD corrections 
which are known |5|, || to modify the shape of the W^-mass peak and the distributions 
of others kinematical variables. Since event shape variables are often used to extract 
WW production from the background it is highly desirable to have a complete next- 
to-leading order (NLO) study of these distributions for WW events. Furthermore, 
calculations of QCD corrections to 0(a 2 a 2 s ) four-jet production have recently appeared 
0. Combining these results with a complete calculation of QCD corrections to all 
C(a 4 ) four-fermion processes it would be possible to obtain NLO predictions for any 
four-jet shape variable at the NLC providing new means of testing perturbative QCD. 

In this contribution we present the complete calculation of QCD corrections to 
e + e~ — > u d ji Up, e + e~ — > u d e u e and e + e~ — > u d s c. . While this is only a first step 
in the calculation of all 0(a 4 ) four-fermion processes at NLO, these reactions include 
the most important source of background, namely single W production, providing a 
natural setting for a first study of the role of QCD corrections to gauge invariant sets 
of four-fermion production diagrams. A gauge invariant description of e + e~ — > u d s c 
requires in the unitary gauge the eleven diagrams shown in fig. 1. This amplitude is 
known as CC11 in the literature. The subset of three diagrams labeled (e) and (f) in 
fig. 1, in which both W 7 s can go on mass-shell, is known as CC03 and is often used for 
quick estimates of WW production. The ten diagrams required for e + e~ — > u d \i u^ 
(CC10) and the twenty diagrams required for e + e~ — > u d e u e (CC20) are similar and 
can be readily deduced from our figures. 

In most instances QCD corrections have been included "naively" with the sub- 
stitution Y w — > T w (l + 2/3a s /ir) and multiplying the hadronic branching ratio by 
(1 + a s /ir). This prescription is exact for CC03 when fully inclusive quantities are 
computed. However it can only be taken as an order of magnitude estimate even for 
CC0 3 in the presence of cuts on jet properties, as discussed in H] 

2 The impact of QCD corrections on the angular distribution of the decay products of a W and 
their application to on-shell VF-pair production is discussed in ref. M 
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It is well known that differential distributions can be more sensitive to higher order 
corrections than total cross-sections in which virtual and real contributions tend to 
cancel to a large degree. It is therefore necessary to include higher order QCD effects 
into the predictions for WW production and decay in a way which allows to impose 
realistic cuts on the structure of the observed events. 




Figure 1: Tree level diagrams for e + e — > u d s c. The dashed lines are W's. 




Figure 2: Basic combinations of loop diagrams. All virtual QCD corrections to electroweak 
four fermion processes can be computed using these three sets. The quark wave functions 
corrections are not included because they vanish, in the massless limit, using dimensional 

regularization. 



2 



Calculation 



One-loop virtual QCD corrections are obtained by dressing all tree diagrams with 
gluon loops. Defining suitable combinations of diagrams as in fig. 2 one can organize 
all contributions in a very modular way, as shown in fig. 3. As already mentioned in 
H all QCD virtual corrections to 0(a 4 ) four-fermion processes can be computed using 
the set of loop diagrams shown in fig. 2. 

The real emission contribution can be obtained attaching a gluon to the quark 
line(s) of the tree diagrams in all possible positions. The required matrix elements 
have been computed using the formalism presented in ref. with the help of a set of 
routines (PHACT) [[KJ which generate the building blocks of the helicity amplitudes 
semi-automatically. 

The calculation of the virtual corrections has been performed in two different ways, 
with identical results. In the first case we have used the standard Passarino-Veltman 
lTj reduction procedure, while in the second we have used the new tecniques presented 



in 
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In order to be able to integrate separately the real and virtual part one has to 
explicitly cancel all singular contributions in each term in a consistent way. To this 
aim we have used the subtraction method. Benefits of this method are twofold. First, 
an exact result is obtained and no approximation needs to be taken; second, all singular 
terms are canceled under the integration sign and not at the end of the calculation, 
leading to better numerical accuracies. This is especially relevant for the present case, 
since we are aiming for high precision results, with errors of the order of a per mil. 
We have found it particularly convenient to implement the dipole formulae of reference 
|I3| . These are a set of completely general factorization expressions which interpolate 



smoothly between the soft eikonal factors and the collinear Altarelli-Parisi kernels in 
a Lorentz covariant way, hence avoiding any problem of double counting in the region 
in which partons are both soft and collinear. 



All integrations have been carried out using the Monte Carlo routine VEGAS P 



In the absence of a calculation of all O(a) corrections to four-fermion processes, we 
have only included the leading logarithmic part of ISR using the (3 prescription in the 
structure functions, where (3 = ln(s/m 2 ) — 1. Beamstrahlung effects have been ignored. 



Results 

In this section we present a number of cross sections and of distributions at -y/i = 500 
GeV. We have used a s = .123 in order to conform to the choice made for the Joint 
ECFA/DESY Study. Initial state radiation is included in all results. 

For the NLC workshop the so called NLC/TH set of cuts have been agreed on: 

• the energy of a jet must be greater than 3 GeV; 

• two jets are resolved if their invariant mass is larger than 10 GeV; 

• jets can be detected if they make an angle of at least 5° with either beam. 
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Figure 3: One loop gluonic corrections to the u d line of the e + e — > u d s c process. 
Similar corrections on the s c line must be included. 



This set of cuts will be referred to as "canonical" in the following. We have pre- 
ferred a different criterion for defining jets which is closer to the actual practice of the 
experimental collaborations. For mass reconstruction studies we have used the Durham 
|T5[ scheme with — 1. x 10" 3 . The four-momenta of the particles which have to be 
recombined have been simply summed. If any surviving jet had an energy smaller than 
3 GeV it was merged with the jet closest in the Durham metric. 



Previous studies [161 have shown that the differences between the total cross sections 



obtained from CC10 and CC11 and those obtained with CC03 are at the per mil level. 
The non resonant background is far more important when there is an electron in the 
final state: the cross section, with canonical cuts, for CC20 at y/s = 190 GeV is larger 
than the cross section calculated from CC03 by several percent. When the full NLO 
results is compared with the naive-QCD (nQCD)predictions, large effects have been 
found in observables like the average shift of the mass reconstructed from the decay 
products from the true W mass for fully hadronic processes. 

In fig. 4 and 5 we show some representative distribution for the semileptonic events. 
In fig. 4a we present the normalized distribution of the angle between the muon and 
the closest jet for CC10. The events in fig. 4a pass all NLC/TH cut with the exception 
of the minimum 9^. With increasing collider energy, the two W's tend to fly apart 
with larger relative momentum and therefore the probability of a jet to overlap with 
the lepton decreases. This behaviour is clearly visible in fig. 4a which shows that jets 
are typically well separated in angle from the charged lepton and that the differences 
between the exact distribution and the one with "naive" corrections is confined to very 
large angles. 

Fig. 4b shows the distribution of the minimum angle between any jet and either 
beam. We have separated the contribution of two-jet and three-jet final states making 



it possible to estimate the effect of different angular cuts on the cross section. 

The total NLO cross section with canonical cuts is .23554(4) pb, to be compared 
with the nQCD results of .23065(3). The two predictions differ by more than 2%. 

Fig. 5a and 5b present the normalized distributions of the angle between the elec- 
tron and the closest jet and of the minimum angle between any jet and either beam, 
respectively, for e + e" — > u d e v e . 

In fig. 6 we compare the NLO spectrum of the average reconstructed W mass for 
CC11 with the nQCD case result. All events with at least four observed jets have 
been retained in the plots. The two candidate masses are obtained forcing all events to 
four jets, merging the two partons which are closest in the Durham scheme, and then 
selecting the two pairs which minimize 



a;, = (m r1 - M w y + (m R2 - M w y . (i) 

where Mm M R2 are the two candidate reconstructed masses and My/ is the input W 
mass. In fig. 6 the dashed line refers to the nQCD results while those of the full NLO 
calculation are given by the continous line. Our plot is obtained using only the basic 
set of cuts described above. 

Fig. 6 shows that at NLO the mass distribution is shifted towards lower masses and 
a long tail for rather small average masses is generated. Because of the larger relative 
momentum of the two W's it is less likely that partons from the decay of one W end up 
close to the decay products of the other VT-boson, therefore the effect is smaller than 
at Lep 2 energies 



In fig. 7 we present the distributions of the following four-jet shape variables \TB 

• the Bengtsson-Zerwas angle: xbz = ^[(Pi x P2), (P3 x P4)] (fig- 7a); 

• the Korner-Schierholz-Willrodt angle : 

§ksw = l/2{Z[(pi x p 4 ), (p a x p 3 )] + Z[(pi x p 3 ), (p 2 x p 4 )]} (fig. 7b); 

• the angle between the two least energetic jets; a 34 = Z[p3,P4] (fig. 7c); 

• the (modified) Nachtmann-Reiter angle: 6* NR = Z[(pi — P2), (P3 — P4)] (fig- 7d). 

The numbering % = 1 . . . , 4 of the jet momenta p^ corresponds to energy-ordered four- 
jet configurations (E\ > E 2 > E 3 > E4). We compare the exact NLO results with the 
distributions obtained in nQCD and with the results obtained at tree level from the 
standard background reactions e + e~ — > q q g g and e + e~ — > q\ q± q 2 q 2 - 

In all subplots of fig. 7 the full NLO results are given by the continous line and the 
nQCD prediction is given by the dashed line. The q q g g and qi q\ q 2 q 2 tree level 
background distributions are given by the chain-dotted and the dotted line respectively. 
The shape variables are computed following the procedure outlined in ref. [p| where the 
Durham cluster algorithm is complemented by the E0 recombination scheme, namely 
if the two particles % and j are merged the pseudo-particle which takes their place 
remains massless, with four-momentum: 

Ei + E- 

E new = Ei + Ej, Pnew = 1 ; r(Pi + Pj)- (2) 

|Pi + Pjl 
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Using this clustering procedure all five-jet events are converted into four-jet events, 
then each event is used in the analysis if minj J=14 ?/jj > y cut with y cut = 0.001. In 
particular no minimum angle between jets and either beams is required. 

It should be stressed that four-jet shape variables in WW events measure the 
correlations between the hadronic decays of the two W's and therefore it should be 
explicitely checked whether existing codes, NLO calculations or parton shower Monte 
Carlo programs, successfully reproduce the experimental curves. 

A separation based on shape-variables of WW events from the background seems, 
at first sight, to be rather difficult at the NLC The variable which most clearly discrim- 
inates between signal and background is the angle between the two least energetic jets 
a 34 . While the signal distribution peaks in the backward direction the background is 
almost flat. On the contrary, the Bengtsson-Zerwas angle and the Korner-Schierholz- 
Willrodt angle distribution from CC11 are almost indistinguishable from those gen- 
erated from the q q g g background. Some additional sensitivity is provided by the 
Nachtmann-Reiter angle distribution. Both signal and background peak at small an- 
gles, but while the former is almost negligible in the backward direction, the latter 
shows a large tail which extends to 180°. The distributions obtained in nQCD are very 
similar to the full NLO results, contrary to what happens at lower energies [[lTj. 

Conclusions 

We have described the complete calculation of QCD radiative corrections to the charged 
current processes e + e~ — > u d \i v m e + e~ — > u d e v e and e + e~ — > u d s c which are 
essential in order to obtain theoretical predictions for 1^-pair production with per mil 
accuracy. The amplitudes we have derived are completely differential, and realistic 
cuts can be imposed on the parton level structure of the observed events. For CC11 we 
have presented the distribution of the average reconstructed W mass in fully hadronic 
events and the distribution of several four-jet shape variables. For CC10 and CC20 the 
distributions obtained in nQCD are very similar to the exact NLO results, however the 
corresponding total cross sections can differ by more than 2%. 
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Fig. 4: Distribution of the angular separation of the \i~ from the closest jet (a) and 
of the minimum angular separation between any jet and either beam (b) at y/s = 500 
GeV for CC10. The continuous, dotted and dot-dashed histograms are exact NLO 
results while the dashed histogram refers to nQCD. 
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Fig. 5: Angular separation of the e~ from the closest jet (a) and between any jet and 
either beam (b) at ^/s = 500 GeV. The continuous, dotted and dot-dashed histograms 
are exact NLO results while the dashed histogram refers to nQCD. 
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Fig. 6: Average mass distribution at y/s = 500 GeV. All NLC/TH cuts are applied. 
The continuous histogram is the exact NLO result while the dashed histogram refers 
to nQCD. 
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Fig. 7: Four-jet shape variables at ^fs = 500 GeV. The full NLO results (continous 
line) is compared with the nQCD prediction (dashed line) and with the tree level 
background distributions from q q g g (chain-dotted line) and q± q~i q2 q<i (dotted line). 




Fig. 7: Four-jet shape variables at ^Js = 500 GeV, continued. 



